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A possible gluon-condensate-induced modified-gravity model with f{R) oc 
|i?|^/^ has been suggested previously. Here, a simplified version is pre- 
sented using the constant fiat-spacetime equilibrium value of the QCD gluon 
condensate and a single pressureless matter component (cold dark matter, 
CDM). The resulting dynamical equations of a spatially-fiat and homoge- 
neous Robertson-Walker universe are solved numerically. This simple em- 
t^^ . pirical model allows, in fact, for a careful treatment of the boundary condi- 

^ ! tions and does not require a further scaling analysis as the original model 

00 \ did. Reliable predictions are obtained for several observable quantities of the 

^ \ homogeneous model universe. In addition, the estimator Eq^ proposed by 

• ' Zhang et al. to search for deviations from standard Einstein gravity, is cal- 

Q ■ culated for linear sub-horizon matter-density perturbations. The QCD-scale 

O . modified-gravity prediction for Eg{z) differs from that of the ACDM model 

by about ±10% depending on the redshift z. 
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1. Introduction 

It has been suggested that the gluon condensate (Shifman, Vainshtein 
& Zakharov, 1979; Narison, 1996) of quantum chromodynamics [QCD] in 
the g-theory framework (Klinkhamer & Volovik, 2008a,b, 2010) may lead to 
a particular type of /(-R) modified-gravity model (Klinkhamer & Volovik, 
2009a). The cosmology of this model was explored in Klinkhamer (2010a). 
Here, the model is simplified even further, allowing for a detailed study of 
the corresponding cosmological equations and eliminating all arbitrariness. 

The model of this article is, however, considered from a different point of 
view than the one used in the original article (Klinkhamer & Volovik, 2009a). 
Here, the aim is not to derive the particular model used from the underlying 
theory (general relativity and QCD), but to consider the model per se and 
to test the hypothesis whether or not this model gives a reasonably accurate 
description of the observed "accelerating universe" (cf. Riess et al., 1998; 
Perlmutter et al., 1999; Komatsu et al., 2009, and references therein). In 
other words, the approach is in the sprit of Kepler rather than Newton. 

Still, the model is not chosen at random from all possible f{R) modified- 
gravity models, but obeys certain physically motivated selection criteria, as 
discussed in Section 2.1. The selected QCD-scale modified-gravity model is 
detailed in Section 2.2, with a few technical points relegated to Section 2.3. 

The corresponding cosmological model with a single matter component 
(cold dark matter, CDM) is introduced in Section 3.1. The resulting dif- 
ferential equations are given in Section 3.2. Three observable quantities of 
the homogeneous cosmological model are then discussed in Section 3.3. A 
fourth observable, this time of the inhomogeneous cosmological model, is 
also considered, namely, the estimator Eq, which was designed (Zhang et 
al., 2007) to characterize the linear growth of matter-density perturbations 
and which has been measured recently for moderate redshifts (Reyes et al., 
2010). In Section 4, the same four observable quantities are calculated ana- 
lytically for the standard ACDM-model universe (Weinberg, 1972; Carroll, 
Press, & Turner, 1992; Sahni & Starobinsky, 2000; Peebles & Ratra, 2003; 
Perivolaropoulos, 2010) and may serve as benchmark results. 

The numerical results for both homogeneous and inhomogeneous observ- 
able quantities of the QCD-scale modified-gravity universe are discussed in 
Section 5 and compared with those of the ACDM universe. Final comments 
are presented in Section 6. 



Even though not necessary for a proper understanding of the present ar- 
ticle, the reader may wish to have a look at Appendix A of a recent re- 
view (Klinkhamer & Volovik, 2011), where the basic idea of g-theory is out- 
lined. 

2. QCD-scale modified-gravity model 

2.1. Selection criteria 

There exists an infinity of f{R) modified-gravity models and an infinity 
thereof reproduces the main characteristics of the observed homogeneous uni- 
verse [see, e.g.. Brans & Dicke (1961); Starobinsky (1980); Bertolami (1986); 
Sotiriou & Faraoni (2010)]. In this paper, a subset of these theories is ob- 
tained by imposing the following conditions: 

1. Empty fiat spacetime is a solution of the field equations resulting from 
the action with the function f{R) replaced by the constant /(O). 

2. The action of the f{R) modified-gravity model involves only known 
energy scales from general relativity and the standard model of ele- 
mentary particle physics [e.g., the energy scales -Eqcd = 0(10^ MeV) 
and -Epianck = ^Jhc^ /{S-k Gn) ~ 2.44 x 10^^ GeV, for the particular 
model of the next subsection]. 

3. The asymptotic de-Sitter solution from the f{R) modified-gravity model 
has only integer powers of the Hubble constant H in the reduced action 
[e.g., the terms (E'pianck)^-^^ and {Eqcb)^H, for the particular model 
of the next subsection]. 

4. The action of the f{R) modified-gravity model has dimensionless cou- 
pling constants roughly of order unity [e.g., the coupling constant rj, 
for the particular model of the next subsection] . 

The first condition is motivated by the need to solve the main cosmological 
constant problem [CCPl] (Weinberg, 1989), namely, why is the energy scale 
of the effective gravitating vacuum energy density pv negligible compared to 
the basic energy scales of the standard model of elementary particle physics 
(not to mention the Planck energy scale). A first step towards solving CCPl 
has been made using g-theory (Klinkhamer & Volovik, 2010). The third 
condition is to have a stationary cosmological solution (de-Sitter universe) 
which is analytic and consistent, as discussed in Section III of Klinkhamer 
& Volovik (2008b). The second and fourth conditions are for simplicity's 



sake, without the need to introduce new physics at ultralow energies [that 
is, energies of the order of meV, as indicated by the observed accelerating 
universe; cf. Riess et al. (1998); Perlmutter et al. (1999); Komatsu et aL 
(2009)]. 2 

Remark that the ACDM model [i.e., the cosmological model from stan- 
dard general relativity with a genuine cosmological constant A and a cold- 
dark-matter (CDM) component] already does not satisfy the first condition. 
Moreover, the required value of A in the ACDM model (Weinberg, 1972; 
Carroll, Press, & Turner, 1992; Sahni & Starobinsky, 2000; Peebles & Ra- 
tra, 2003; Perivolaropoulos, 2010) is completely ad hoc (the meV scale men- 
tioned above) and leaves CCPl hanging in the wind, which is theoretically 
unsatisfactory, as explained in, e.g., the second paragraph of Section I of 
Mukohyama (2004). Just to avoid any misunderstanding, the flat ACDM 
model (Weinberg, 1972; Carroll, Press, & Turner, 1992; Sahni & Starobinsky, 
2000; Peebles & Ratra, 2003; Perivolaropoulos, 2010) is perfectly satisfactory 
experimentally but not theoretically. The QCD-scale modified-gravity model 
proposed in the next subsection has the potential to be incorporated in a 
theoretically satisfactory framework (g-theory) but, first, needs to be shown 
experimentally satisfactory. 

2.2. Representative model 

The following simplified modified-gravity action (Klinkhamer & Volovik, 
2009a; Klinkhamer, 2010a) satisfies the criteria of Section 2.1 and will be the 
starting point of the present article: 

5grav,o = [ d'x^/^^\KoR{x)-r]\R{x)\'/'{qof^ 

+Cm[^{x)]], (2.1) 

with h = c = 1 from the use of natural units, the gravitational coupling 
constant Kq = (IGttGo)"^ > 0, the dimensionless coupling constant r] > 



^Of all four conditions, the third has admittedly the weakest theoretical motivation. 
But additional arguments may perhaps come from future results on the quantum-field- 
theoretic de-Sitter state, which may hold some surprises in store [see, e.g., Polyakov (2010) 
and references therein!. 



[standard general relativity has 77 = ] , the constant equilibrium gluon con- 
densate go = (-E'qcd)^, and the generalized matter field ip{x) which includes 
the fields of QCD (gluons and quarks) and further possible fields responsi- 
ble for the observed CDM component of the present universe."^ Throughout, 
the conventions of Weinberg (1972) are used such as the metric signature 
( — h ++), except for the Riemann tensor, which is taken to have a further 
minus sign. Specifically, the Ricci curvature scalar of a de-Sitter universe is 
given by -R = 12 H^ > with Hubble constant H. 

It remains to be seen if there exist other f{R) modified-gravity mod- 
els which satisfy the criteria of Section 2.1 and describe the observed uni- 
verse equally successfully. Perhaps a model exists involving the electroweak 
scale (ArkaniHamed et al., 2000; Klinkhamer & Volovik, 2009b; Klinkhamer, 
2010b, 2011). Returning to the energy scale of QCD, remark that the second 
term of the integrand in (2.1) corresponds to the linear term {EqqdYH men- 
tioned in the second criterium of Section 2.1 and that this linear term may be 
considered to be the leading term of a power series in H (higher-order terms 
will be discussed shortly). Anyway, model (2.1) based on the QCD energy 
scale is a perfect example (essentially unique up to higher-order terms) of a 
model satisfying the criteria of Section 2.1. As explained in Section 1, the 
aim of the present article is to study the use of this model as a compact 
description of the observed universe. 

In the standard formulation of modified-gravity models where the Ricci 
scalar R of the Einstein-Hilbert action density is replaced by i? + f{R), the 
proposed modification is given by 

f{R) = -\R\'/yLo, (2.2a) 

with length scale 

Lo = r/-iJfo(go)~'/'. (2.2b) 

It is important to state, right from the start, that the theory (2.1) is only 
considered to be relevant over cosmological length scales (small curvatures 



■^As discussed in Klinkhamer & Volovik (2009a), the vacuum energy density of the 
QCD fields with dynamics governed by Cm in (2.1) is compensated by an appropriate 
value go of the gluon-condensate q(x) field, resulting in a vanishing effective cosmological 
constant. Other (non-QCD) contributions to the gravitating vacuum energy density are 
perhaps canceled by the self-adjustment of similar g-type fields or by a different mechanism 
ahogether; cf. Ftn. 2 of Khnkhamer & Volovik (2009a). 



\R\)- For smaller length scales (larger curvatures |-R|), other terms than the 
single root term (2.2) may become important. One example would be having 
f{R) of (2.2) replaced by the following extended function: 

/ext(i?) = -^_^^|^|l/2°^^, (2.3) 

with another positive coupling constant ( of order unity. This particular 
function differs from the one given by Klinkhamer (2010a) in that its second 
derivative is positive for all finite values of R. 

2.3. Additional remarks 

The model of the previous subsection with a constant value go of the 
gluon condensate is a simplified version of the model with a dynamic gluon 
condensate field q{x) considered in Klinkhamer (2010a). It was found in 
(Klinkhamer, 2010a) that the gravitating vacuum energy density pv{(l) oc 
(? ~ QoY of fhs corresponding cosmological model is negligible during the 
late evolution of the model universe and that the dynamic vacuum variable 
q{x) is effectively frozen to its constant equilibrium value go with pv{(lo) = 0. 
For this reason, it makes sense to restrict the consideration to model (2.1) 
with just a constant value go- 

The actual value of the gluon-condensate go entering the nonanalytic 
modified-gravity term (2.2) arises from the flat-spacetime part of the La- 
grange density. Cm, and is taken to have a standard (positive) value go ~ 
(300 MeV)^; see Shifman, Vainshtein & Zakharov (1979); Narison (1996). 
Later, it will be shown that this particular value for go implies 77 ^ 2 x 10~^. 
Observe that, strictly speaking, the quantity go of the second term in the 
integrand of (2.1) need not be equal to the equilibrium value of the gluon 
condensate g(x), but can, in principle, be given by any typical QCD energy 
density. Still, the identification of go as the equilibrium gluon condensate 
is maintained throughout this article, because, with this identification, the 
coupling constant 77 in (2.1) is defined unambiguously. 

The relation between the gravitational constant Gq = G{qo) = {16ttKo)~^ 
from (2.1) and Newton's constant On (Cavendish, 1798; Mohr, Taylor & 
Newell, 2008) is rather subtle, but, in this article, the approximate equality 
Go ~ Gn is simply taken to hold (see the next paragraph for the argument). 
For these numerical values of the constants 77 (go)^'^^ and Go, the length scale 



entering the modified-gravity action term (2.2) has the value Lq ~ 10^® 
3 Gpc. 

The field equations from (2.1) are fourth order and it is advantageous to 
switch to the scalar-tensor formulation, which has field equations of second 
order but involves an extra scalar field (Brans & Dicke, 1961; Starobinsky, 
1980; Sotiriou & Faraoni, 2010).^ The scalar-tensor theory equivalent to 
(2.1) has been given by Eq. (2.2) of Klinkhamer (2010a), with q replaced 
by go and changing the sign in front of the (j) R term, where (p G (— c>o, 1) 
is the dimensionless Brans-Dicke scalar field. The corresponding potential 
t/(0) oc rj'^/{l — 0) is such that it leads to the chameleon effect (Khoury 
& Weltman, 2004) and the effect is taken to be operative for Cavendish- 
type experiments on Earth (Cavendish, 1798; Mohr, Taylor & Newell, 2008), 
giving Go ~ Gjv (see also Endnote [39] in Klinkhamer (2010a) for further 
discussion) . 

3. One-component QCD-scale modified-gravity universe 

3.1. Pressureless matter component 

In order to study cosmological aspects of the modified-gravity model (2.1), 
consider a spatially fiat {k = 0) Robertson-Walker metric (Weinberg, 1972) 
with scale factor a^r) for cosmic time r (later, t will denote the dimensionless 
quantity) and matter described by homogeneous perfect fiuids. Recall that 
the Hubble parameter is given by H{t) = [da{T)/dT]/a{T). 

In fact, only one matter component (labeled n = 2) will be considered in 
the rest of this article, namely, a perfect fiuid of pressureless nonrelativis- 
tic matter [e.g., cold dark matter] with energy density Pm2{t) and constant 
equation-of-state (EOS) parameter wm2 = Pkn/ Pm2 = 0. This is a simpli- 
fication of the model considered by Klinkhamer (2010a), which had also a 
dynamic gluon-condensate component (labeled n = 0) with EOS parameter 
Wq = —1 and an ultrarelativistic-matter component (labeled n = 1) with 



'^ The modified Einstein equation [see, e.g., Eq. (6) of Sotiriou & Faraoni (2010)] also 
has singular terms proportional to ji?]"^/^ or higher negative powers of \R\. It remains to 
be seen whether or not this leads to unacceptable behavior. As mentioned in Klinkhamer 
(2010a), appropriate versions of the model may satisfy solar-system tests because of the 
chameleon effect (to be discussed shortly). See also the related discussion in the penulti- 
mate paragraph of Section 5. 



EOS parameter wmi = 1/3- In order to allow for an easy comparison with 
the results of Klinkhamer (2010a), the label n = 2 will be kept for the single 
matter component considered in this article. 

Apart from this cold-dark- matter component, the only other input of the 
cosmological model is the modified-gravity term from (2.1). Using the scalar- 
tensor formalism, there is then the auxiliary Brans-Dicke scalar 0(r) to con- 
sider, without direct kinetic term but with the nontrivial potential U{(f)) 
already mentioned in Section 2.3. 

3.2. Dimensionless variables and ODEs 

The following dimensionless variables t, /i, /, m, s, and r can be introduced: 



r E 


= tK,/{r^{q,fl'), 


(3.1a) 


H{r) E 


E h{t)r^{q,f"/K,, 


(3.1b) 


t/(r) E 


= u{t) r^^q^f' / Kl , 


(3.1c) 


0(r) E 


E s(t). 


(3.1d) 


pA/2(r) E 


= rM2{t)v\q^f'VK^^ 


(3.1e) 



where, different from Klinkhamer (2010a), the rescaling is done with the com- 
bination If] {q^Y^^/KQ = 7] {qoY^^ (IGvrGo) appearing in the modified-gravity 
term (2.2a). All dimensionless quantities are denoted by lower-case Latin 
letters. 

From the action (2.1) in the scalar-tensor formulation (Klinkhamer, 2010a), 
there is the following closed system of 4 first-order ordinary differential equa- 
tions (ODEs) for the 4 dimensionless variables h{t), s(t), v(t), and rM2it): 

h = -2/,2-i^, (3.2a) 



6 ds 



(3.2b) 



1 , 2 1 an , , 

V = -rM2-^hv--u + -s—-, (3.2c) 

6 3 3 as 

rM2 = -3/irM2, (3. 2d) 

where the overdot stands for differentiation with respect to the dimensionless 



cosmic time t and and the dimensionless Brans-Dicke potential u is given by 

With the solution of the above ODEs for appropriate boundary conditions, 
it is possible to verify a posteriori the Friedmann-type equation, 

h^s + hv={rM2-u)/6, (3.4) 

which, in general, is guaranteed to hold by the contracted Bianchi identities 
and energy-momentum conservation (Weinberg, 1972). Hence, if the solution 
of ODEs (3.2) satisfies (3.4) at one particular time, then (3.4) is satisfied at 
all times. 

The boundary conditions for ODEs (3.2) are obtained by setting t = tstart 
in the following functions: 

.^PP'--(t) = (l-^t + ^A , (3.5b) 

^approx^^^ = s'^PP™^(t), (3.5c) 

-l?r (t) = lUl-'-^t). (3.5d) 



These functions provide, in fact, an approximate solution of the ODEs (3.2) 
and generalized Friedmann equation (3.4) in the limit t — ;■ 0. 

Purely mathematically, remark that the ODEs (3.2) have no free param- 
eters and that the boundary conditions are fixed completely by (3.5) for 
sufficiently small tstart- Physically, however, tgtart must be small enough but 
still larger than the time tcross corresponding to the QCD crossover at a 
temperature Tcross ~ Eqcb ~ 300 MeV [see Klinkhamer (2010a) for further 
discussion]. Specifically, one must take a tgtart value obeying 1 ^ tstart 3> 

tcross ~ V -EQCoZ-Epianck ~ 10"^^ for f] ~ 10"^. 



3.3. Observables quantities 

In order to test the cosmology resulting from QCD-scale modified-gravity 
model (2.1), four observable quantities will be considered (Klinkhamer, 2010a; 
Zhang et al., 2007): tph(tp), Wx(tp), -Zmfioctl^j, tp), and Eq"^"{z), where z 
stands for the redshift with respect to the present epoch t = tp to be defined 
below. As a preliminary, note that the generalized Friedmann equation (3.4) 
gives 

uJx + uJM2 = 1 , (3.6a) 

UJX = rx/(6/j^) =CJy +Cjgrav, (3.6b) 

uJv = rv/{Qh^) = 0, (3.6c) 

cJgrav = 1 — s — s/h — u/{6 h?) , (3.6d) 

wm2 = rM2/{Qh'^)- (3.6e) 

Without genuine vacuum energy density [rv{q) = KqT]''^ (Qo)'^^"^ Pv{q) = 
for q = qo], the only new ingredient in (3.6) is cJgrav, as it vanishes for the 
standard theory with u = and s = 1. Remark that rx of (3.6b) does not 
correspond to a real physical energy density but is a mathematical quantity 
inferred from the variables uy and cJgrav defined by (3.6c) and (3.6d). 

The energy density parameters of (3.6) are defined in terms of the grav- 
itational constant Gq = G{qo), which may, in principle, be different from 
Newton's constant G^ as measured by laboratory experiments on Earth (Ca- 
vendish, 1798; Mohr, Taylor & Newell, 2008). For this reason, the parameters 
have been denoted by u and not by the standard symbol Q. Still, the pre- 
cise relation between Go and Gn will be important only once in Section 5, 
namely, when the absolute age of the universe is calculated and, then, the 
equality Go = Gn is taken to hold [assuming the chameleon effect (Khoury 
& Weltman, 2004) to be relevant for earth-based experiments, as discussed 
in Section 2.3]. Now, turn to the four observables mentioned above. 

The first observable is simply the expansion rate of the universe in units 
of the inverse of its age, 

tph{tp)=td{t)/a{t)l^^^, (3.7) 

where a{t) is the scale factor of the Robert son- Walker metric. This observ- 
able quantity has been evaluated for the "present epoch," which is taken to 



10 



be defined by the moment t = tp when uJM2iip) = 1/4 or Pxijp)/ PAnijp) = 3 
from (3.6a). The fiducial value ujM2{tp) = 1/4 will be used in the following 
and is more or less consistent with the available data [see, e.g., Komatsu et 
al. (2009)]. 

The second observable corresponds to the effective EOS parameter of the 
unknown component X, whose model value can be extracted from (3.2) and 
(3.4): 

_ , , 2 f aa 1 

Wx[tp) = 



3 \{a)'^ 2 J 1 - ujm2 

u + Ahs + 2s 

u + 6hs — rj\/2 (1 — s) 



t — t-n 



(S.t 



Again, this observable quantity has been evaluated at the present epoch t = 
tp. For later times, the right-hand side of (3.8) shows that wx of the modified- 
gravity model (2.1) approaches the value —1 in the limit of vanishing matter 
content rM2 and constant Brans-Dicke scalar s as t — > oo. 

The third type of observable follows from the transition of deceleration to 
acceleration. In mathematical terms, this time corresponds to the nonsta- 
tionary infiection point of the function a{t), that is, the value tj at which the 
second derivative of a{t) vanishes but not the first derivative. Specifically, the 
infiection point t = ti corresponds to the following redshift for an observer 

diXi L ^ Zp -^ Li '. 

2;mflect(ii, tp) = a{tp)/a{ti) - 1 . (3.9) 

In order to prepare for the fourth observable of the empirical model (2.1), 
turn to the linear growth of sub-horizon matter-density perturbations in the 
Newtonian gauge (Song, Hu & Sawicki, 2007). For small enough matter- 
density-perturbation amplitudes Am2 = 5fM2/fM2 and large enough wave- 
lengths (but still within the horizon), the following dimensionless linear ODE 
needs to be solved: 

AM2(t) + 2 hit) AM2{t) - 7T ^77 rM2{t) AM2{t) = , (3.10) 

3 s[t) 

which is equivalent to, for example, Eq. (13) of Tsujikawa et al. (2009) or 
Eq. (15) of Lombriser et al. (2010), with an extra factor (4/3) (1/0) multiply- 
ing the original non-derivative term. See also de la Cruz-Dombriz, Dobado 
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& Maroto (2008), where the approximate ODE (3.10) has been tested for an 
|i?|^/^ modified-gravity model. 

From the combined solutions for a{t) and AM2{t), the following linear 
growth parameter is obtained by differentiation: 

/3^^^ = 1#^. (3.11) 

a In a h Am2 

For the early phase with (/) = s ~ 1, /i ~ 2/3 1~^, and rM2 ~ 6 h!^, the linear 
growth parameter calculated from (3.10) is 



/3early= (V33-l)/4^ 1.186. (3.12) 

The corresponding matter-density perturbation is given by 

which provides an approximate solution of the ODE (3.10) in the limit t — ?■ 0. 
Now, consider the estimator Eq introduced by Zhang et al. (2007) with 
the purpose of searching for deviations from Einstein gravity. According to 
the original paper of Zhang et al. (2007) and also the paper of Lombriser et 
al. (2010), the theoretical expression for Eq{z) in the modified-gravity model 
(2.1) is 

with /ij = df /dR. The last expression of (3.14) is given in terms of the scalar- 
tensor formalism, where the relevant field equations (Klinkhamer, 2010a) 
have been used. From the numerical solution of the ODEs to be presented 
in Section 5, the values of (j){z) and (3{z) are readily obtained and, thereby, 
the numerical value of Eq^°°{z). In the early phase with ~ 1 and /5 given 
by (3.12), one has 

^g''"L»1 ~ ^M2(tp)//3early ^ 0.2108 , (3.15) 

for the fiducial value uJM2iip) = 0.25 mentioned earlier. 

4. Benchmark results from the ACDM model 

Four observables of the QCD-scale modified-gravity model (2.1) have been 
discussed in Section 3.3 and will be evaluated numerically in Section 5. In 
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this section, the corresponding values are given for the spatially flat ACDM- 
model universe (Weinberg, 1972; Carroll, Press, & Turner, 1992; Sahni & 
Starobinsky, 2000; Peebles & Ratra, 2003; Perivolaropoulos, 2010). 

Recall that the relevant spatially fiat ACDM-model is completely defined 
by the value of the cosmological constant A and the condition that flM2 
equals 1/4 at present. Here, the standard matter energy- density parameter 
in terms of Newton's constant G^ is given by Weinberg (1972) 

^M2{Tp) = PM2{rp)/Pcrit{rp) , (4.1a) 

Pcrit(rp) = 3H\Tp)/{8nGN), (4.1b) 

where the present epoch occurs at cosmic time t = Tp. As mentioned in 
Section 2.1, the ACDM model may be theoretically unsatisfactory but has 
been found to give an excellent description of the observed accelerating uni- 
verse (Riess et al., 1998; Perlmutter et al., 1999; Komatsu et al., 2009). For 
this reason, it can provide benchmark results to compare other models with. 
Analytic results for the first three observables of Section 3.3 are as fol- 
lows (Klinkhamer, 2010a): 



ACDM 



\TpH{Tp), Wx{Tp), ^inflect(Ti, Tp 

= ((4 arcsinh 73)7(3 73), -1, (6^/^ - 1)) 

'l.Ol, -1, 0.817) . (4.2) 

The fourth observable is simply (Zhang et al., 2007) 



^t'i^) LcDM = ^M2{rp)/PACDM{Z) , (4.3) 

where /3acdm(^) can be calculated analytically from the exact solution of 
the linear ODE corresponding to (3.10); see, in particular, the C2 term in 
Eq. (6.67) of Mukhanov (2005). The estimator (4.3) approaches the constant 
value fiM2{Tp) for increasing redshift z, as (3\cdm{z:) goes to 1 for 2; — )• 00. 

5. Numerical results for the modified-gravity universe 

Figure 1 displays the numerical solution of the first-order nonlinear ODEs 
(3.2) with boundary conditions (3.5) at t = tstart = 10~^. The simultaneous 
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numerical solution of the second-order linear ODE (3.10) has been obtained 
for boundary conditions from the approximate solution (3.13) at t = tstart- 
Specifically, (3.13) provides the initial derivative Aj\/2(tstart) for a given ini- 
tial value of AA./2(tstart)- Numerical solutions of (3.10) with other values for 
the initial derivative Ajv/2(^start) have been seen to rapidly approach the ap- 
proximate solution (3.13), provided the dimensionless cosmic time t remains 
small enough. 

The linear t scale of Fig. 1 is convenient for the late evolution of the model 
universe, because, as will be seen shortly, t corresponds to the cosmic time 
r measured in units of approximately 10^° yr = 10 Gyr. A logarithmic scale 
is, however, more appropriate for the early phase and Fig. 2 shows that the 
numerical results from the QCD-scale modified-gravity model (2.1) reproduce 
the Friedmann-Robertson-Walker-type expansion with Hubble parameter 
h = 2/3t~^ and linear growth parameter /3 given by (3.12). Table 1 displays, 
moreover, a stable behavior of the t = tp observables within the numerical 
accuracy (at the one-per-mill level or better, for the quantities shown). 

The results for the first three observable quantities of Section 3.3 are as 
follows: 



TpH{Tp), WxiTp), ZinRcctiTi, Tp) 

p^ (0.917, -0.662, 0.523) . (5.1) 

These results may be compared to the ACDM-model values (4.2). 

The obtained values (5.1) are consistent with those of Table I in Klinkhamer 
(2010a) for the modified-gravity model with a dynamic q field, but with- 
out the need to consider the limit of the mathematical parameter Z = 
(go)^''^ Kq^ = 2 (-EQcoZ-Spianck)^ to a numerical value of order 10~^^, as Z 
has been scaled away completely in the present simplified model. The only 
place where this hierarchy parameter Z enters is for the dimensional age of 
the present universe and its expansion rate. 

In fact, with Go = G^ and the elementary-particle-physics result (Shif- 
man, Vainshtein & Zakharov, 1979; Narison, 1996) for the flat-spacetime 
gluon condensate go = (-E'qcd)^ = (300 MeV)^, the values {tp, hp) = (1.374, 
0.6673) from the numerical solution (Table 1) give the following results for 
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the present age and expansion rate of the model universe: 

^ 13.2 Gyr, (5.2a) 

Hp = hp7]K^\qof' 

= hpri{l67iGN){EQCDf 

^ 68.1kms~^ Mpc-\ (5.2b) 

r] ^ 2.40xl0-\ (5.2c) 

where the numerical values given in (5.2a) and (5.2b) use the t] value given 
in (5.2c). The modified-gravity model (2.1) can thus give reasonable values 
for both the present age and expansion rate, which traces back to the fact 
that the numerical value for the product tp hp has been found to be close 
to 1. Only the actual value of r] depends on the precise determination of 
go and Gq/Gn- All other results of this paper do not depend on r] , qo, or 
Go as these quantities can be scaled out according to (3.1). The rest of this 
section focusses on two such observable quantities, H{z)/H{0) and Eg{z), 
both considered as a function of the redshift z. 

The behavior of the Hubble parameter H{z) relative to the calculated 
present value (5.2b) is given by the second column of Table 2. The corre- 
sponding modified-gravity results for H{z) agree with the observations re- 
ported in Fig. 9 and Table 2 of Stern et al. (2010), even though the scatter of 
the data is still substantial. The third column of Table 2 gives, for compar- 
ison, the Hubble-parameter ratios from the ACDM model, which also agree 
with the current H{z) observations (Stern et al., 2010) if the measured value 
of H{0) is used as input. 

The numerical solution of Fig. 1 at t = 1.0011 gives the redshift z ~ 0.32, 
the scalar-field value = s ~ 0.7666, and the linear growth parameter 
P ^ 0.7459. These numbers combined with the fiducial value UM2p = 0.25 
result in the following numerical value for the gravity estimator (3.14): 

^g'^1.=o.32^ 0-437, (5.3) 

which agrees within one sigma with the experimental result E^^ = 0.392 ± 
0.065 from Reyes et al. (2010) for a sample of galaxies with an average redshift 
{z) = 0.32. 
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The same quantity E^*^" at 2; ~ 0.3 has also been calculated by Lombriser 
et al. (2010) for an f{R) modified-gravity model with 1 + //j = ~ 1 and 
Bq parameter (Song, Hu & Sawicki, 2007) of order 0.2, giving a value around 
0.35. In order to compare to this result, the value of the Bp parameter (in 
the notation with the present epoch at r = Tp) has been calculated for the 
QCD-scale modified-gravity model (2.1): 



R Irr ( 1 dR 



fn \R din a 
1 



rd\na\ 
V dH ) 


R H\ 




rh) 


T = Tp 



Bp 



2 2Lo^/\R\- 1 
^ 0.246, (5.4) 

where the Ricci scalar is given by i? = 6 {dH/dr + 2 H"^) and the function 
f{R) is defined by (2.2). With this Bp parameter and simply omitting the 
factor 1/(1 + /ij) or 1/0 in (3.14), our result for Eg at z = 0.32 would be 
approximately 0.34, which would agree with the result of Fig. 4 in Lombriser 
et al. (2010). 

Nevertheless, the factor 1/(1 + //?) is unarguably present in the Eq expres- 
sion (3.14) and the correct estimate from the QCD-modified-gravity model is 
(5.3), which is larger than that of Lombriser et al. (2010) but still consistent 
with the direct measurement (Reyes et al., 2010). Observe that the value 
/r(^p) ~ 0.7259 — 1 f« —0.2741 is consistent with the cosmic-microwave- 
background data according to Table III of Lombriser et al. (2010), but ap- 
parently inconsistent with the cluster-abundance data according to Table IV 
of the same reference. However, as mentioned before, the QCD-modified- 
gravity model is assumed to hold only for the very largest length scales and 
not galaxy-cluster length scales (or, a forteriori, solar-system length scales). 
Interestingly, the cluster-abundance data can be used to constrain the ex- 
tensions of the simple |-R|^/^ term, one possibility having been mentioned in 
(2.3). 

Finally, additional values for Eq^'^" at selected redshifts are given in the 
fourth column of Table 2. This table also compares the QCD-modified- 
gravity results for E^'^° with those of the ACDM-model universe. [The 
ACDM value for Eq^'^° ai z = 0.32 is 0.396, whereas the modified-gravity 
value has already been given in (5.3) and the experimental value just below 
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that equation.] Putting the theoretical predictions for Eg{z) in Table 2 
next to the simulated data in Fig. 1 of Zhang et al. (2007) suggests that it 
may be difficult for future surveys (e.g., the Square Kilometer Array on the 
ground or the Joint Dark Energy Mission and the Euclid Satellite in space) to 
distinguish between the two theoretical models but perhaps not impossible. 

6. Discussion 

In this article, a simple empirical model has been proposed with a QCD- 
scale modified-gravity term (2.1) and a single pressureless matter component 
(cold dark matter, CDM). This particular f{R) modified-gravity model has 
been selected on physical grounds (Section 2.1), but is, in the end, solely 
used as an efficient way to describe the main aspects of the late evolution of 
the universe, having only two fundamental energy scales, -Eqcd and -Epianck, 
and a single dimensionless coupling constant, rj. 

With the elementary-particle-physics value for the equilibrium gluon con- 
densate go = {Eqcy^Y = (300 MeV)^, the measured age of the universe 
fixes the dimensionless coupling constant t] of model (2.1) to the value of 
approximately 2.4 x 10~^. As emphasized in Klinkhamer (2010a), the effec- 
tive coupling constant t] may ultimately be calculated from QCD and general 
relativity; see also Schiitzhold (2002); Bjorken (2004); Urban & Zhitnitsky 
(2010a,b); Holdom (2011) for further discussion of the possible relation of 
QCD and dark energy. The connection of the QCD-scale modified-gravity 
model with the g-theory approach to solving the main cosmological constant 
problem has already been mentioned in Section 2.1. 

The go and r] values quoted in the previous paragraph result in a model 
postdiction for the Hubble constant (5.2b), which is within 10% of the ob- 
served value. Theoretically, the last mathematical expression for the Hubble 
constant in the middle of (5.2b) is quite remarkable, as it involves only New- 
ton's constant Gn and the cube of the energy scale -Eqcd- Independent of 
the go and rj values, there are further model predictions for the dimension- 
less quantities (5.1) and (5.3), which are again in the same ball park as the 
observed values. The modified-gravity model prediction of —0.7 for the quan- 
tity Wx{Tp) as defined by (3.8) [compared to the ACDM value of —1] may 
perhaps provide a crucial test, as long as independent measurements of the 
present values of H, dH/dr, and pcdm = Pm2 can be obtained. 

Moreover, model predictions for the redshift dependence of the gravity 
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estimator Eq have been given in Table 2. At the very largest length scales 
(comoving wavelengths of the order of hundred Mpc or more but still less 
than the horizon scale) and relatively low redshifts z ~ 0.75, the QCD- 
modified-gravity results for Eq differ by some +10% from the ACDM-model 
values. The QCD-modified-gravity results for Eg differ by some —10% from 
the ACDM-model values for z ^ 3, but it is not yet clear how these redshifts 
can be probed observationally. 

As it stands, the QCD-scale modified-gravity model (2.1) gives a remark- 
able description of the main aspects of the late evolution of the universe. 
With many forthcoming ground-based experiments and space missions, fu- 
ture measurements oiWx{Tp) and Eg{z) may suggest alternative f{R) func- 
tions, such as the one given by (2.3). Of course, if experiment indeed finds 
evidence for nonzero f{R), it is up to theory to provide the proper under- 
standing. 
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Figure 1: Numerical solution of the modified-gravity cosmological ODEs (3.2), with Brans- 
Dicke scalar potential (3.3) and pressureless matter having a dimensionless energy density 
tmi [see (3.1) for the definitions of the dimensionless variables used]. The boundary con- 
ditions follow from the approximate solution (3.5) evaluated at tstart = 10^^. The initial 
value a(istart) for the scale factor a(t) is taken as 10""^ but has no direct physical rele- 
vance. The linear matter-density-perturbation ODE (3.10) is solved simultaneously with 
boundary conditions from (3.13) evaluated at tgtart- The figure panels are organized as 
follows: the panels of the first column from the left concern the scale factor a(<) and 
the Hubble parameter h = {da/dt)/a, those of the second column the Brans-Dicke scalar 
s{t) [the QCD-modified-gravity model (2.1) being studied in the scalar-tensor formalism], 
those of the third column the matter energy density rM2 [the bottom panel of this column 
showing the linear growth parameter /? of subhorizon matter-density perturbations] , and 
those of the fourth column derived quantities [the bottom panel showing the gravity esti- 
mator E^'^° defined by (3.14)]. The three panels of the fourth column can be combined 
to give the behavior of h{z)/h{0) and Eq{z) shown in Table 2. The several energy-density 
parameters TD and the effective "dark-energy" equation-of-state parameter wx are defined 
in (3.6) and (3.8), respectively. 
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Figure 2: Semi- log plot of numerical quantities from Fig. 1, with /3cariy defined by (3.12). 
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Tabic 1: Function values for the "present epoch" [defined by uJM2{tp) = 0.25] from the 
numerical solution of the modified-gravity cosmological ODEs (3.2), with Brans-Dicke 
scalar potential (3.3) and boundary conditions (3.5) taken at different values of igtart- The 
values for the linear growth parameter /3 of subhorizon matter-density perturbations follow 
from the simultaneous numerical solution of ODE (3.10) with boundary conditions from 
(3.13) evaluated at tgtart- 
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Table 2: Numerical modified-gravity results from Fig. 1 for the Hubble parameter H{z) 
relative to its present value -ff(O) and the gravity estimator Ea{z) defined by (3.14), 
taking uJM2{tp) ~ 0.25 for the Eg{z) values quoted. The Eq value for redshift z = 10^ 
is essentially equal to the analytic result (3.15). For comparison, also values arc given for 
the spatially flat ACDM model with ^M2{Tp) = 0.25 and Eg{z) deflned by (4.3). 
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